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Abstract 

We determine the gluonic field configuration sourced by a heavy quark undergoing 
arbitrary motion in M = 4 super- Yang-Mills at strong coupling and large number of 
colors. More specifically, we compute the expectation value of the operator Tr[F 2 +. . .] 
in the presence of such a quark, by means of the AdS / CFT correspondence. Our re- 
sults for this observable show that signals propagate without temporal broadening, 
just as was found for the expectation value of the energy density in recent work by 
Hatta et al. We attempt to shed some additional light on the origin of this feature, 
and propose a different interpretation for its physical significance. As an application 
of our general results, we examine (Tr[F 2 + ...]} when the quark undergoes oscilla- 
tory motion, uniform circular motion, and uniform acceleration. Via the AdS/CFT 
correspondence, all of our results are pertinent to any conformal field theory in 3 + 1 
dimensions with a dual gravity formulation. 
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1 Introduction and Summary 
1.1 Motivation 

When a charge moves, it produces a propagating disturbance in the associated gauge 
field. The problem of determining the spacetime profile of this disturbance for an 
arbitrary charge trajectory was solved long ago for classical electrodynamics [I], but, 
under various guises, remains of significant interest today in the context of quantum 
non-Abelian gauge theories. In recent years, gauge/gravity duality [21 OH] has given 
us a useful handle on this and many other problems for a varied class of gauge theories 
in the previously inaccessible regime of strong coupling, via a drastic and surprising 
rewriting in terms of string-theoretic degrees of freedom living on a curved higher- 
dimensional geometry. 

Among the known examples of the duality, the best understood subclass is that of 
conformal field theories (CFTs), where the relevant curved geometry is asymptotically 
anti-de Sitter (AdS). In this paper we will use this AdS/CFT correspondence to study 
the propagation of disturbances in the gluonic field produced by a moving heavy quark 
in a strongly-coupled conformal gauge theory, in the limit where the number of colors 
is large. For concreteness, we will phrase our analysis in terms of M = 4 super- Yang- 
Mills (SYM), even though the results we will obtain are equally relevant to other 
CFTs in 3 + 1 dimensions, can easily be extended to CFTs in other dimensions, and 
might also be expected to apply at a qualitative level in some of the non-conformal 
examples of the gauge/gravity correspondence. 

The AdS/CFT correspondence states that a heavy quark moving in the vacuum 
of M = 4 SYM is dual to a string moving on a pure AdSs geometry. More pre- 
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cisely, the quark corresponds to the endpoint of a string, whose body codifies the 
profile of the non-Abelian (near and radiation) fields sourced by the quark. When 
the quark/endpoint moves, it generally produces a wave running along the body of 
the string, which corresponds to the wave in the gluonic field whose profile we are 
interested in determining. 

The translation between string and gauge theory disturbances was first explored 
in [5], which used tools developed in (31 E] to study the dilaton waves given off by 
fluctuations on an otherwise static, radial string in AdSs, and infer from them the 
profile of the dual gluonic field observable (Tr F 2 (x)) in the presence of an oscillating 
quark. Under the assumption that these oscillations are small, the authors of [5] 
treated the string dynamics in a linearized approximation. Their results painted an 
interesting picture of wave propagation in M = 4 SYM: in contrast with the standard 
Lienard-Wiechert story of the classical Abelian case, where signals propagate strictly 
at the speed of light, the waves in (Tr F 2 (x)) were found to display significant temporal 
broadening, just as one would expect given that points on the non-Abelian field 
arbitrarily far from the quark can themselves reradiate. 



Figure 1: Schematic representation of emission of dilaton waves by a string in AdS. The 
radial direction of AdS is denoted by z and runs vertically, with the AdS boundary (z = 0) 
located at the bottom of the figure. Any given point x on the boundary receives a dilaton 
wave from each point along the string. As a result, the gauge theory wave is a superposition 
of components with all possible time delays. 

As shown in Fig. 1, this feature emerges naturally in the gravity side of the 
correspondence, because motion of the string endpoint generates waves that move up 
along the body of the string, and each point on the string then emits a dilaton wave 
that travels back down to the observation point x on the AdS boundary, where, via the 
AdS/CFT recipe for correlation functions [3], the value of (Ti F 2 (x)) is deduced by 
assembling together all such contributions. The dilaton waves originating from points 
on the string that are further away from the AdS boundary give rise to components 
of the SYM wave that have a larger time delay. This makes perfect physical sense, 
because, through the UV/IR connection [7], such points are known to be dual to 
regions of the gluonic field further away from the quark. So, even when treating the 
string waves at the linearized level, the very fact that we are dealing with a string 
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leads to gauge theory disturbances that display the nonlinear propagation expected 
from the non-Abelian character of the strongly-coupled SYM fields. 

The simple observable (TrF 2 (x)) that was the focus of [5] served well to exhibit 
the main features of propagating disturbances in the gluonic field, but did not allow 
a definite identification of waves with the characteristic l/|x| 2 falloff associated with 
radiation, i.e., contributions that transport energy to infinity. Of course, in the case 
of an accelerating quark, fluctuations in the near-fields are fully expected to be ac- 
companied by radiation proper. The unambiguous detection of the latter calls for 
examination of the SYM energy-momentum tensor, (T Mi/ (x)), which in the gravity 
side of the correspondence requires a determination of the gravitational waves emit- 
ted by the string. Indeed, some time after [5], the l/|x| 2 falloff was established in the 
same setup through a calculation of the (time- averaged) energy density, (T 00 (x)) [8]. 

In more recent years, motivated by potential contact with the phenomenology 
of the quark-gluon plasma [9], analyses of both (TiF 2 (x)) and (T^x)) have been 
carried out in the case of a heavy quark moving at constant velocity through a thermal 
plasma (where an exact solution for the corresponding string embedding is available 
[HJj [TT]). in a large body of work that includes [12] and has been reviewed in [13j . 
The results are again compatible with the expected nonlinear dynamics of the (in this 
case, finite-temperature) SYM medium. 

Given these antecedents, it came as a surprise when, back at zero temperature, 
additional calculations going beyond the linearized string approximation found, first 
for special cases [HI [15] and then for an arbitrary quark trajectory [16], that the 
ensuing energy density (Tqo(x)) displays no temporal broadening, and is in fact as 
sharply localized in spacetime as the corresponding classical profile. As shown in 
[TBI US], the fact that disturbances in (Too (a;)) propagate strictly at the speed of 
light does not actually conflict with the UV/IR argument presented three paragraphs 
above, because the full gluonic profile can be understood as arising purely from a 
contribution of the string endpoint, and this is the reason why it features a single 
time delay. 

Our work was motivated by the tension between these two sets of results for the 
gluonic fields in vacuum. How can it be that the (TrF 2 (x)) profile obtained in [5] 
displays significant temporal broadening, while the (T o(x)) pattern deduced in [IS] 
does not? A natural strategy to further explore and attempt to resolve this tension 
is to carry out the (Tt F 2 (x)) calculation beyond the linearized string approximation, 
and for an arbitrary quark trajectory, to put it on a par with the computation in |16j . 
This is what we set out to do in this paper. 

1.2 Outline and main results 

We begin in Section[2]by presenting the ingredients and recipe for our calculation. The 
desired one-point function of TrF 2 (x) (or, more precisely, of the Lagrangian density 
operator ([S])), in the presence of a moving quark, must be extracted via ([7]) from the 
leading near-boundary behavior of the dilaton field (TTS]) sourced by the string profile 
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(lT6"j) . The latter is the unique string embedding codifying the gluonic fields generated 
by a quark undergoing arbitrary motion, under the assumption that such fields are 
retarded or purely outgoing, i.e., they propagate outwards from the quark to infinity. 
This solution was obtained in [17] for an infinitely heavy quark and generalized in 
p~8j [T9"l 120] to the case of a quark that has a finite mass, and, consequently, a finite 
size (i.e., Compton wavelength) z m given by ([2]). 

The actual computation is carried out in Section [3J A key feature is that, when 
expressed in terms of appropriately geometric variables (the quark proper time r 
and the invariant AdS distance Q), a total derivative is found to appear in the 
worldsheet integrand (1231) , allowing one of the two integrals to be done trivially. This 
feature results from a nontrivial cancelation (noted below (|2~T|) ) that in turn originates 
from the retarded structure of the string embedding ( 1T-41) . An analogous cancelation 
was found in the energy density computation in [IS], arising there from an interplay 
between the various components of the string energy-momentum tensor. Since our 
analysis involves only scalar quantities, it becomes clear that the cancelations in 
question are not intrinsically tied to the tensorial nature of the (Tqq(x)) derivation. 

Our final result for the gluonic profile (TrF 2 (x)) in the case of a quark with 
finite mass is somewhat involved: it is written in (T5TT) in terms of auxiliary (tilde) 
variables whose explicit dependence on the quark's velocity, acceleration, jerk and 
snap, as well as on the external force (and the first and second derivatives thereof) 
it is subjected to, is given in (132]) - (I34]) . (l36]) - (l37j) . It is easy to check that, when the 
quark is static, our general expression correctly reproduces the known [21] finite-mass 
result (14*21) . For arbitrary motion, the gluonic field at any given observation point is 
found to depend only on dynamical data evaluated at a single retarded time along 
the quark trajectory, specified by (1381) or, in non-covariant form, ([39]) . This is the 
same surprising characteristic demonstrated for the radiation component of the field 
in [151 [16]. 

In Section HI we apply the general result (13T1) to three specific examples. The 
first is the oscillating quark that was the focus of [5] and the main motivation for 
our work. We start by performing a numerical analysis of the case (not covered by 
[5]) where the quark has a finite mass, and display the resulting gluonic profile in 
Fig. 4. Moving on to the case of infinite mass, we show explicitly that the integrated 
expression obtained in [5], which manifestly displays the gluonic field as a sum of 
contributions with all possible time delays, actually gives the same result (within its 
range of validity) as the linearized and infinitely-massive version of our final formula 
(T3T[) for the gluonic profile, which incorporates a single time delay. This proves that, 
despite appearances, there is in fact no conflict between [5] and [151 US]: it is just that 
the choice of worldsheet coordinates and the lack of an exact solution for the string 
embedding prevented the authors of [5] from being able to carry out the integral 
explicitly to find the net retardation pattern. 

Our second example is uniform circular motion, studied in Section 14.21 to com- 
plement the energy density analysis of [13]. Our third and final example, examined 
in Section I4.3[ is uniform acceleration, for which [221 1231 [23] had already made some 
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interesting physical inferences based on the structure of the string embedding. 

In our final Section we go back to our result for arbitrary quark motion, and 
discuss its physical implications. We consider first the case of an infinitely heavy 
quark, which was the only one contemplated by the (TrF 2 (x)) computation in [5] 
and the (Too (a;)) calculations in [TU EOS [16] . in this limit, the quark becomes pointlike 
(z m — > 0), and our result simplifies drastically, taking the form ( H3j) . For an arbitrary 
spacetime trajectory of the quark, this is manifestly just the boosted Coulomb profile 
associated with a uniformly translating quark with the retarded position and velocity 
inferred by projecting back along the past lightcone of the observation point. The full 
profile agrees (up to an overall constant) with the corresponding Lienard-Wiechert 
result in classical electrodynamics [I]. The observable (TrF 2 (x)) determined in this 
paper is thus seen to yield an image of the gluonic profile for a pointlike quark that 
is in complete consonance with the one obtained via (T 00 (x)) in [TB] , with the main 
difference being that, just like in electrodynamics, the former is sensitive only to 
disturbances in the near field of the quark, whereas the latter also explicitly incor- 
porates radiation. Both here and in [16] signals in the gluonic field are found to 
propagate strictly at the speed of light (as indicated below (JI2D), with no temporal 
(or, equivalently, radial) broadening. 

Our more general result (13TI) . relevant for a quark with finite mass, differs from 
([43]) and from the energy density obtained in [16] in two ways. First, the retarded 
quark data are read off at a source point that is related to the observation point 
through the timelike interval (l38j) . In other words, signals are found to propagate at 
a (variable) subluminal speed, which as seen in (1631) can in fact be arbitrarily small. 
Second, the field profile depends on more data than just the position and velocity of 
the quark. As explained in Section [5], both of these features are naturally associated 
with the fact that the quark is no longer pointlike. 

Even with these differences, it is still true that in our general result fl3T|) the field at 
a given observation point is controlled by data at a single retarded event, as in [T5l[l~6"]. 
The discrepancy between this pattern and the one reported in [5] was discussed in 
[T5] in terms of a fundamental distinction between the quark's near field (which is 
all that is visible in {TrF 2 (x))) and its radiation field (which dominates (Tqq(x)) at 
long distances), but our results show that there is in fact no such distinction, because 
both components display exactly the same unbroadened propagation. The reason on 
the gravity side is the same in both cases: the fact that (in accord with the UV/IR 
connection) the result can be expressed in terms of a contribution arising purely 
from the string endpoint. Moreover, we know from Section H~T1 that the discrepancy 
between [5] and [151 EE] is om y apparent, and that, in spite of the fact that the gluonic 
field emerges as a superposition of contributions with all possible time delays, the net 
result evidences only the smallest of these delays. 

In [15l [16] it was argued that the lack of broadening is unphysical and points 
to a deficiency of the AdS/CFT calculation. The argument visualizes the radiation 
process in terms of emission of gluons that are themselves able to reemit, and would 
be expected to do so profusely at strong coupling. This would bring into play large 
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quantities of off-shell gluons, which would naturally be expected to propagate at 
subluminal speeds and therefore to lead inevitably to temporal/radial broadening 
of the emitted field. This 'parton branching' picture is essentially a perturbative 
rephrasing of the discussion in [5] (recalled above) of non-Abelian reemission by the 
gluonic field at all different length scales, and has been shown to be consistent with 
several other AdS/CFT results [25J. Given the apparent incompatibility of this picture 
with the lack of broadening, the authors of [TSJ US] went on to suggest that, contrary 
to widespread belief, the supergravity approximation to physics on the gravity side 
does not capture the full quantum dynamics of the large A and large N c limit, and 
tried to identify the missing element as arising from longitudinal fluctuations in a 
heuristic lightcone gauge calculation. 

As explained in Section [51 even though we find in this paper that (TrF 2 (x)) dis- 
plays exactly the same retardation pattern as (T 00 (x)), we do not subscribe to the 
point of view of [T5| [T6] . As seen in the calculation for arbitrary quark trajectory in 
Section [3] (as well as in [16]), and also when we make contact in Section H~T1 between 
[5] and our general result ( l3Tj) . the supergravity description assembles the gluonic field 
precisely in the physically expected manner, by summing over contributions reradi- 
ated from all possible length scales (as depicted schematically in Fig. 1), associated 
with all possible time delays. It is therefore not true that the absence of net temporal 
broadening in the final AdS / CFT result implies that the supergravity approximation 
is somehow leaving out the expected non-Abelian rescattering. Rather, we interpret 
the no-broadening result as a prediction of the AdS/CFT correspondence for the net 
pattern of propagation in the CFT at strong coupling and with a large number of 
colors. 

We stress in particular that the appearance in ( 123]) of a total derivative in the 
worldsheet integrand, which enables us to present the resulting dilaton field as a 
pure endpoint contribution, does not mean that points on the body of the string 
do not contribute, but only that their aggregated contribution can be reexpressed 
in terms of the behavior at the edge of the integration region. The same applies 
then on the SYM side: the fact that the gluonic field at the observation point can 
be reported purely in terms of the behavior of the quark at a single retarded time 
does not indicate that only that instant contributes, but only that the cumulative 
effect of summing over the contributions from all previous emission events (arising 
via non-Abelian rescattering at all possible length scales) can be reexpressed in terms 
of the aforementioned behavior. This is explicitly shown by our calculations, by 
those of [TU [151 dSj, and also by a recent surprising reformulation of the latter as 
a superposition of gravitational shock waves emitted by each point along the string 

It is important to keep in mind that the form of our result depends crucially on the 
retarded structure of the worldsheet embedding ( f!6|) . which in turn follows from the 
assumption of a purely outgoing condition for the gluonic field generated by the quark. 
For any choice other than the purely outgoing (or purely ingoing) SYM configuration, 
we would expect not to obtain a total derivative on the worldsheet, and this would 
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then lead to a final gluonic profile showing finite temporal/radial broadening. At least 
in retrospect, it is natural for the purely outgoing condition to impose a restriction 
on the overall retardation pattern of the field, because arbitrary reradiation from all 
points on the non-Abelian medium would result in wave scattering back towards the 
quark. That the final gluonic profile at any given observation point, which receives 
contributions from the entire quark trajectory, can be reexpressed in terms of data at 
a single retarded time, is also not surprising per se. At the calculational level, that 
happens every time we evaluate an integral in terms of the behavior at the integration 
endpoint. What is remarkable about the restriction predicted by AdS/CFT is that 
it involves only a finite number of quark data at the relevant instant. 

In short, then, we believe our results resolve the apparent conflict between [5] and 
[T5l fT6] , and support an interpretation of the lack of temporal broadening that does 
not challenge the ability of the supergravity approximation of AdS/CFT to capture 
the full dynamics in the strong-coupling limit. 



2 Ingredients of the Computation 

We are interested in studying the gluonic field sourced by a heavy quark in a strongly- 
coupled gauge theory. Gauge/gravity duality grants us access to many different se- 
tups, and in particular, to conformal field theories (CFTs) in any dimension, but 
for concreteness we will focus on M = 4 super- Yang-Mills (SYM) with gauge group 
SU(N C ). This is a conformally invariant theory with a vector field, 6 real scalars and 
4 Weyl fermions, all in the adjoint representation of the gauge group. The AdS/CFT 
correspondence [2] asserts that this theory, on (3 + l)-dimensional Minkowski space- 
time, is fully equivalent to Type IIB string theory on the Poincare patch of the 
AdS 5 x S 5 geometry, 

ds 2 = G mn dx m dx n = — (-dt 2 + dx 2 + dz 2 ) + R 2 dnl , (1) 

z 1 

(with a constant dilaton and N c units of Ramond-Ramond five-form flux through 
the five-sphere). The radius of curvature R is related to the SYM 't Hooft coupling 
A = 9y M N c through 



where l s denotes the string length. In more detail, the state of IIB string theory 
described by ([T|) corresponds to the (symmetry-preserving) vacuum of the gauge the- 
ory, and the closed string sector describing small or large fluctuations on top of 
it encodes the gluonic (+ adjoint scalar and fermionic) physics. The coordinates 
x^ = (t, x) parallel to the AdS boundary z = are directly identified with the gauge 
theory spacetime coordinates, the radial direction z is mapped to a variable length 
(or, equivalently, inverse energy) scale in SYM [TJ, and the five-sphere coordinates 
are associated with the global SU(A) internal (R-) symmetry of SYM. These an- 
gular coordinates will play no role in our analysis, so our results will hold equally 
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well in the more general case where the S 5 is replaced by a different compact five- 
dimensional space X 5 , which corresponds to replacing M = 4 SYM with a different 
(3 + l)-dimensional CFT. 

The introduction of an open string sector associated with a stack of Nf D7-branes 
in the geometry <^ is equivalent (27] to the addition, on the gauge theory side, of 
Nf hypermultiplets (each composed of a Dirac fermion and 2 complex scalars) in 
the fundamental representation of the SU(N C ) gauge group, that we will refer to as 
'quarks'. For Nf <C N c , we are allowed to neglect the backreaction of the D7-branes 
on the geometry]]] in the gauge theory this corresponds to working in a 'quenched' 
approximation that ignores quark loops. The D7-branes cover the four gauge theory 
directions x^, and extend along the radial AdS direction up from the boundary at 
z = to a position where they 'end' (meaning that the S 3 C S 5 that they are wrapped 
on shrinks down to zero size), whose location z = z m is related to the mass m of the 
quarks through 

«- = • (2) 

An isolated quark is dual to an open string that extends radially from the location 
z = z m on the D7-branes to the horizon of the Poincare patch, z — > oo. We will 
describe the dynamics of our string in first-quantized language, and since we take 
it to be heavy, we are allowed to treat it semiclassically. In gauge theory language, 
then, we are coupling a first-quantized quark to the gluonic (+ other SYM) field(s), 
and carrying out the full path integral over the strongly-coupled field(s) (the result 
of which is codified by the AdS spacetime), but treating the path integral over the 
quark trajectory x(t) in a saddle-point approximation. 

In the nonperturbative framework provided to us by the AdS/CFT correspon- 
dence, a quark with finite mass (z m > 0) is automatically not 'bare' but 'composite' 
or 'dressed'. This can be inferred, for instance, from the expectation value of the glu- 
onic field surrounding a static quark, as in [21] and the calculations we will perform 
in the next section, or from the deformed nature of the quark's dispersion relation 
[j~8| [T9| [20] . The characteristic thickness of the 'gluonic cloud' surrounding the quark 
is given by z m , which is thus the analog of the Compton wavelength for our non- 
Abelian source. 

Most of the time we will write the quark trajectory in Lorentz-covariant notation 
parametrized by proper time, x(r) = x^(r). In our analysis below we will have need 
to refer to the velocity, acceleration, jerk and snap of the quark, which will be denoted 
as 

dx _ dv -> da _ dj 

v = Tt' a =dt' J = ^' s = ^' (3) 

or, in four-vector form, 

dx dv da dj 

v = — , a = — , j = — , s = —, (4) 
ar dr dr dr 



1 More precisely, the relevant condition is XNf <g; N c [28] 
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which by definition satisfy the relations 



v =—1, v-a — 0, v ■ j = —a . (5) 

By the AdS/CFT dictionary, it is the endpoint of the string on the D7-branes that 
directly corresponds to the quark, whereas the body of the string codifies the profile 
of the gluonic field sourced by the quark, which is precisely the information we are 
after. Specifically, our aim is to probe the gluonic field configuration by determining 
the one-point function of the operator 

F 2 = — g— Tr {F^F^ + $j][$ 7 , $ J ] + fermions} , (6) 

where & , / = 1, . . . , 6, denote the scalar fields of the Af = 4 SYM theory. Op* is 
essentially the SYM Lagrangian density. As a shorthand, from this point on we will 
refer to it (including the normalization constant in front) as simply TrF 2 . Since, as 
mentioned above, we work in the approximation where the quark follows a definite 
trajectory, its presence amounts to the addition of the corresponding Wilson line, due 
to which we expect the one-point function of Tr F 2 to be nonvanishing. The latter 
operator is known to be dual to the IIB (s-wave) dilaton field 4>{x,z) [221 ED], and 
the standard GKPW recipe for correlation functions [3] at large N c and A relates 
its one-point function to a variation of the supergravity action with respect to the 
boundary value of 0. The connection can be summarized as [5] 

(TrF 2 (x)) = -\im(^dMx,z)) , (7) 



in terms of a rescaled dilaton field ip = Vt^R^ifi^K 2 . 

Our task, then, is to calculate the dilaton field sourced by the string dual to the 
quark, and to pick out the 0(z A ) term in its expansion near the boundary of AdSs. 
In the linearized approximation appropriate at large N c , <p is obtained simply by 
convolving the string source with the retarded dilaton propagator [B] 



D(U) 



d 



4ir 2 R 3 dU LVl - U 2 



2U 2 - 1 



e(i-|z/|) 



which depends only on the invariant AdS distance between the observation (un- 
primed) and source (primed) point, 

{t-t'f- (x-x) 2 - (z-z'f _ (x - x'f + z 2 + z' 2 

2z~z~< ~ 2z~I> ' {) 

The symbol inside the brackets of (jHJ) denotes the Heaviside step function, which 
implements causality, as we will elaborate on below. 

The string dynamics is prescribed as usual by the Nambu-Goto action 



Sng - 

In 



^ J d 2 a v 7 - det g ab (10) 
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where g a b = d a X m dbX n G mn (X) (a,b = 0,1) denotes the induced metric on the 
worldsheet. We work for the time being in the static gauge a = t, a 1 = z, so the 
string embedding is described by X(t,z). Given that the string endpoint represents 
the quark, the latter's trajectory is read off from the string profile via 

x(t) =X(t,z m ) . (11) 

Since we have mentioned that for z m > our fundamental color source is not pointlike, 
we should be more precise: the AdS/CFT dictionary identifies X(t, z m ) as the location 
of the most UV contribution to the gluonic field, and it is this that we adopt as a 
natural definition for the position of the quark. This definition has proven to be useful 
to capture the physics of a range of phenomena including Brownian motion (31] and 
(as we will recall below) radiation damping [191 ED], but, given the extended nature 
of the quark, is of course not unique (see, e.g., [32]). 

Using (jSj), and taking into account the coupling (in the Einstein frame) between 
the dilaton and the string worldsheet, one finds the relation [5] 



f d 



2U 2 — 1 

■B(1-\U\) 



(12) 



Together with (JTj), this will yield the result we are seeking. To put this to use, we are 
missing only one ingredient: we must know the explicit form of the string embedding. 

Of course, just like the specification of a quark trajectory does not uniquely de- 
termine a gluonic field configuration, knowing the motion of the string endpoint is 
not enough to select a unique string profile. Additional information is needed, in the 
form of initial or boundary conditions. For a quark in vacuum, the configuration of 
most evident interest is the retarded one, where waves in the gluonic field move out 
from the quark to infinity, rather than the other way around, or some (nonlinear) su- 
perposition thereof. Luckily, the corresponding solution of the Nambu-Goto equation 
is known, for an arbitrary time-like quark/endpoint trajectory. In the case where the 
quark is infinitely massive (z m = 0), this solution can be written as [17] 

X(t r ,z) = x(t r )+ 7==^= , (13) 

y/i - v(t r y 

t ytf ; Z ^ — t>p I 



V 71 - v(t r 

where the behavior of the string at a given time t and radial depth z is seen to be 
completely determined by the behavior of the quark/string endpoint at the earlier, 
retarded time t r H The definition of t r implicit in (fT3"j) can be shown to follow from 
projecting back to the AdS boundary along a curve that is null on the string world- 
sheet [TTJ, in analogy with the Lienard-Wiechert story in classical electrodynamics. 
In Lorentz-covariant language, the solution (fT3l takes the simple form 

X»(t, z) = x"(r) + zv"{t) , (14) 



2 Reversing the sign of the z-dependent terms in (|14j) yields instead an advanced solution |17j . The 
form of the most general Nambu-Goto solution for an arbitrary endpoint trajectory is not known. 



10 



where r is the quark proper time corresponding to t r jf| From now on we will mostly 
work with this parametrization of the string embedding, where the natural notion of 
proper time r associated (modulo a rescaling) with the endpoint has been extended 
to the full worldsheet by following the upward null geodesies. We should stress that 
r is associated with t r and is therefore still a retarded time parameter, in spite of the 
fact that, for brevity, we are not labeling it with the subindex r. 

For the case of a quark with finite mass (and therefore, finite size), z m > 0, the 
string endpoint is at z = z m , and is subject to the boundary condition ffTTj) . As in 
[T8| [T9l [20] , the associated string embeddings can be regarded as the z > z m portions 
of the solutions (fl4|) . which are parametrized by data at the AdS boundary z = 0. 
From this point on we will use tildes to label these (now merely auxiliary) data, 
and distinguish them from the actual physical quantities (velocity, proper time, etc.) 
associated with the endpoint /quark at z — z m , which will be denoted without tildes. 
In this notation, (Tl4]) reads 

X"(f,z) =x*(?) + zv fl (f) . (15) 

As shown in [191 120], this can be rewritten purely in terms of physical (z = z m ) data 
as 

A-(r, z) = x»{t) + {Z ~ ^^J^ , (16) 

where = (2it / y/XjJF^ with JF^ = (— 7-F • v,jF) the external four-force that needs 
to be exerted on the quark to get it to follow the prescribed trajectory^ Expression 
f fTEj) clearly reduces to ( 1T3]) as z m — > 0. As explained in [20], the inclusion of the 
force at finite z m is a convenient way of summarizing a dependence of the string 
embedding on an infinite number of higher derivatives of the quark trajectory, whose 
appearance is natural, given that the dressed quark is an extended object. The 
AdS/CFT correspondence makes it possible to deduce that the trajectory and the 
force are connected through the equation of motion 




^ 2 ( d i^-Sr^ 

W [ 1 - 



T» ~ -KZ-T 1 T 1 , (17) 



which incorporates the effects of radiation damping and is thus a nonlinear general- 
ization of the classic Lorentz-Dirac equation, codifying (inside the parentheses on the 
left-hand side) a non-standard dispersion relation for the dressed quark, as well as 
(in the second term on the right-hand side) a Lorentz-covariant formula for its rate 
of radiation PUH]]. 



3 To avoid confusion, we emphasize here that r is defined with the Minkowski metric appropriate 
for the SYM theory, and therefore differs from the proper time of the string endpoint (computed 
with the AdS metric (jTJ) ) by a factor of R/z m (where we are currently discussing the case z m = 0). 

4 On the gravity side of the correspondence, this is achieved by subjecting the string endpoint to 
an electromagnetic field on the D7-brancs, and = —F vll d T x v is then the usual Lorentz four-force. 
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In the following section, we will use the string embedding ( !T6|) in (j!2p to determine 
the resulting dilaton profile, and then extract from (j7j) the desired expectation value 
of the gluonic field generated by the quark. 



3 Gluonic Profile for Arbitrary Quark Motion 

We now proceed to the calculation. If we employ the more physical parametrization 
(r, z) of the string worldsheet, with r the (retarded) quark proper time, the string 
profile takes the form (TT6]) . Computing the induced metric on the worldsheet, one 
finds that 

y/-g(r',z') = . R (18) 

We have found it more convenient to carry out our entire computation with the world- 
sheet parametrized by (f , z), with f the proper time of the auxiliary (and fictitious) 
string endpoint at z — 0, under which the embedding takes the simpler form (TT5]) . 
implying 

V-9(r',z') = ^ (19) 

(which naturally coincides with the z m — > limit of f fT8|) ). and similarly, all subsequent 
expressions turn out to be more compact. Once we obtain a final result for the 
expectation value of Tr F 2 , we will show how to eliminate all auxiliary (tilde) variables 
in favor of the physical (nontilde) ones. 

Given ( JT9i) and (j2J), to determine the dilaton field ( 1T21) we need to compute the 
double integral 



, \ vx r ,~, r dz' d 

tp[x, zj = / get 



f' 

16tt 3 ./-oo " z' 2 dU 



2U 2 — 1 

-Mi-\u\) 



VT^^u 2 



(20) 



where the invariant distance (Q with the source point taken on the embedding (|T5|) 
adopts the form 

LT _ (x - X{f, z')) 2 + z 2 + z' 2 _ (x- x(f')) 2 + z 2 (x- s(f ')) ■ v(f) 
2zz' ~ 2zz' z ' 

Notice how the z' 2 term of the numerator, present in the definition ([9]), has exactly 
canceled due to the z'-dependence of the retarded solution (115)) . This cancelation will 
play a crucial role in the form of our results, and we will return to it in Section |5j 

The presence of the [/-derivative in ( l20l) suggests a change of variables z' — > U . 
From ( 12T1) we see that 

dU = -dM (X - i( ^ + Z * \ , (22) 
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so we are left witb@ 

8tt 3 

where we have defined 

^min 



dr' z 



(x — x(r')) 2 + z 2 

df'z 
(x — x{f')) 2 + z 2 

(x - x(f')) 2 + z 2 



dU 



d 



f.-, 



2U 2 



2U 2 



^[Y^Tj 2 



dU Lx/T 

e(i-|c/|) 



u 2 

^max 



8(1 -M) 



{x - 5(f)) • 6(f) 



-5(f)) 



(23) 



(24) 



am 



It is now necessary to understand the physical import of the 0(1 — factor. 
Its presence in (|S]) implies that the propagator D{U) has the geometric structure 
shown in Fig. [2j a given source point (x' , z') induces a nonvanishing dilaton on and 
inside its forward lightcone, up to the event where this lightcone reflects back from 
the boundary. Beyond that, there are some observation points (x, z) that lie in 
the future of the source point but are not influenced by it. The region U > 1 
corresponds to observation points that are in the causal past of the event (x' , z') or 
are spacelike related to it. Causality dictates that these points will be beyond the 
region of influence of the source. Exclusion of points at U < — 1, on the other hand, 
is a priori unexpected, because they do obey the naive causality restriction, and are 
only disallowed due to conditions at the boundary of AdS. These are points that can 
be reached by a timelike curve originating at (f , z'), but not by a timelike geodesic. 

In Fig. Owe show the way in which our new coordinates (f , U) cover the string 
worldsheet, which we had originally parametrized with (t',z'). For concreteness, the 
coordinate grid has been plotted for the case where the quark is static; for more 
general motions the overall pattern is similar but somewhat distorted. As mentioned 
above, the constant f (or equivalently, constant t' r ) curves are null; on the (t 1 , z') plane 
they are simply straight lines with slope I/7. The range of constant-f/ curves that 
intersect a given constant-f line runs from U = U m - m {f'\ x, z), where the intersection 
takes place at t' = t' r , z' — 0, to U = f/ max (f ; x, z), which is asymptotic to the line 
and only intersects it at t', z' — > 00. From this it becomes clear that, for any finite f 
and any finite observation point, the point U = U max lies outside of the causal swath 
—1<U < 1 allowed by the Heaviside function. Consequently, only the U = U m \ n 
endpoint contributes to ff23l) . leaving us with 



<p(x,z) 



df' z 



(x — 5(f)) 2 + z 2 



2U 2 



1 



2 

rnin 



e(i-\u n 



(25) 



5 Upon carrying out the U integral to be left with the surface term, one might worry about the 
fact that the quotient in front of the Heaviside function diverges at the endpoints of the interval 
that the latter defines. This concern can be dispelled by noting that 9(1 — \U\) = 0(1 — U)Q(l + U) 
appears within the U derivative, and upon differentiation via the Leibniz rule, would give rise to 
delta functions that precisely cancel these divergences. 
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Figure 2: Penrose diagram of AdS (with boundary directions other than the time 
omitted), showing a Poincare coordinate grid. The shaded yellow region \U\ < 1 
indicates the observation points (x, z) that can be influenced by the given source 
point (x 1 , z'). 



To carry out the remaining integral, it is convenient to perform a final change of 
variables f — > U m - m - From fl24l) we can deduce, via implicit differentiation, that 



zz dU ■ 

T ~~z m + (x- x(f)) ■ (v(f) + z m a(f)) ' 

but proceeding further requires explicitly inverting (|24p . which is impossible to do 
exactly for an arbitrary trajectory. We must remember from (j7J), however, that we 
are ultimately only interested in the value of ip(x, z) in the limit z — > 0, up to order 
z 4 . In this limit the causal swath —1<U < 1 depicted in Fig. [3] in fact collapses 
to a curve on the (t', z') plane, but the function inside the brackets in (125]) still 
varies wildly over the allowed U m i n range, making it necessary to really carry out the 
integral instead of just evaluating the integrand at any point within the interval and 
multiplying times the interval width. A natural strategy is then to Taylor-expand all 
expressions in powers of z. Using ff26l) in ( 125]) . we see that we already have a factor of 
z 2 up front, so we only need to carry out the expansions to order z 2 , because terms of 
order z 3 and higher will not contribute to the one-point function in the z — > limit. 
We start by proposing that 

f = f + hz + f 2 z 2 + 0{z 3 ) , (27) 
where fj = fj(x, U m - m ) are independent of z. Using (}2"7]) we can deduce that 

z(f) = x(f ) + hv(f )z + (f 2 v(f ) + \fl~a{t»)) z 2 + 0(z 3 ) , (28) 
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Figure 3: For a static quark, this is the way in which our new coordinates (f', U) 
cover the string worldsheet which we had originally parametrized with (t',z'). A 
detailed explanation is given in the paragraph before ( l25|) . For ease of visualization, 
we have chosen the observation point at x^ = 0, z = 2 (even though we are ultimately 
interested in z — > 0). The causal swath \U\ < 1 is the region between the 2 dotted 
(red) hyperbolae on the left (which collapses to a line as z — > 0). 



and similarly for v(f'). With this information in hand, we can Taylor-expand ff2^|) in 
a power series in z, match terms of the same order and solve recursively to find the 
coefficients fj. At leading order we get a condition determining fo, 

(x - x(f )) 2 = 2z m (x - x(f )) • v(f ) . (29) 

Proceeding to higher order one finds subsequently 

U ■ z 

1 z m + (x - x') ■ (v' + a'z m ) 

1 Ul in zi(l + (x-x l )-(a' + ~fz m )) 

r 2 



2 [z m + (x- x') ■ (v 1 + a'z m )} 2 [z m + (x - x') ■ (& + a'z m )\ 



3 ' 



where from now on it is understood that functions denoted with a prime are evaluated 
at the retarded time f = f defined implicitly by ( f2"9~T) . 

Using fr2"6]) -( l3"0]) in f l25|) . one can Taylor expand in powers of z, and carry out the 
integral over U m i n to determine <p(x, z). The master formula ([7]) then leads to the 
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gluonic profile we were after, 



(TtF 2 (x)} 



32tt 2 (( x - x>) . v'f [z m + {x- x') ■ {v' + d'z m )f 
x {2 ((2 - x') ■ v'f + z 3 m (l + (x- x') ■ d'f 

+z m ({x - x') ■ v'f [2 + (x - x') ■ (2d' - 0z m - ~s'z 2 m ) - d' 2 z 2 m ] 
+z 2 a (x -x')-v' 4 + 4 ((x - x') ■ d'f + 2(x - x') ■ j'z m 

1 o // ~/\ ~/\ 2 2 I ~!\ ~l 2 ~/2 2 

+3 {(x - x ) ■ j ) z m - (x - x) ■ s z m - a z m 



(31) 



+(x - x') ■ a' (8 + (x - x') ■ (2f - s'z m )z n 



The calculation is straightforward but a bit messy, so we defer the details to the 
Appendix. 

Now that we have our desired result, all that remains is to rewrite the auxiliary 
variables x, v, d, j and s, which describe the motion of a fictitious endpoint at 
z = 0, in terms of the real physical variables x, v, a, j and s, associated with the 
quark/endpoint at z = z m . This connection has been worked out in [20]. Equation 
(25) of that paper states that d translates into nontilde variables according to 



a" 



- 4^ 2 



(32) 



where T 11 is the rescaled version of the external four-force applied on the quark that 
we introduced below ( TTEj) . From this and equation (20) of [2U] it also follows that 



- z 2 T» 



J 72 



(33) 



which in turn implies, via ( TT5|) . 



x M = x 



Comparing (1181) and ( fl9l) we see that 

df = 



1 Z mJ~ 2 



dT 



(34) 



- 4^ 2 

which enables us to deduce from fl32l) that 



(35) 



(i - ztf*) 



(1"4^ 2 ) 



if (36) 
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and 



l - ztT 2 P l + 



2: 



^2 



((T-T)+T)(l-z^)P 



1 - 4^)3/2 



(J" •#•)(! 



4^ 2 ) 3/2 



2™ 



2^: 



l-^OF-J-) 



-4^1-^/. (37) 



where dots over the four-force denote r derivatives. 

Employing ([32}, ([33]), flM]), (SB} and ([37]), we can express the profile fl3TJ) of the 
gluonic field directly in terms of the physical (nontilde) variables that describe the 
quark trajectory. This is our main result. We refrain from writing the resulting 
expression, because it is long and not particularly enlightening. 

All dynamical quantities in the final result are understood to be evaluated at the 
retarded proper time To defined by translating ( 129]) into physical variables, i.e., 



(x - x(r )) 2 = -z. 



(38) 



This equation describes a two-sheeted hyperboloid about the observation point x, 
which is intersected by the quark worldline twice, once on each sheet. By causality, 
the root of interest is of course the one in the sheet to the past of x, which, in 
noncovariant notation, corresponds to the retarded time 



t 



ret 



t 



a/ \x - x(t vct )) 2 + z 



(39) 



This past hyperboloid describes all the points x' = (t ret ,x(t ret )) that can influence 
the given observation point. We can of course read (1381) the other way around, as a 
statement that the events x that can be influenced by a given source point x' = x(tq) 
on the quark worldline, lie on the future hyperboloid at constant timelike interval z m 
from x'. 

Before moving on to applications and a discussion on the physical content of our 
result, let us perform a check on it, by examining what it implies for a free quark. 
In the absence of external forcing, the equation of motion (TTT1) naturally implies that 
the quark moves at constant velocity (although, interestingly, the converse is not true 
[20]). From (!3"2l - (l3~Tj) . we have x = x — z m v, v = v,a=j = s = 0, and therefore (!3T!) 
reduces to 



(TrF 2 (x)) 



V\ 2[(x- x') ■ vf - 4 [(z - x') ■ vf z m + 6 [(x — x') ■ v'} z 



3z: 



32tt 2 



((x — x') ■ v' — z m ) \{x — x') ■ vf 



(40) 

By Poincare invariance, it suffices to evaluate this in the frame where the quark is at 
rest at the origin. Using ( 139]) . we then have (x — x') ■ v' = — (t — t TCt ) = — a/x 2 + z^, 
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so we are left with 



y/X Ax 2 z m + 7z^ + (2x 2 + %z 2 m ) 



(TrF*(x)) = ^— 2 m m ' ' m ^ (41) 

This can be further massaged into the form 

(42) 

which precisely coincides with the gluonic profile determined in [21] for a static quark 
with finite mass. 

For future use we also note that in the limit z m — > 0, where the quark becomes 
pointlike (as well as infinitely massive), our general result ( 1ST]) simplifies drastically, 
and we are left with the single term 



(TrF 2 (x)> = ^5- j , (43) 

16^ [(x-*(to))- W (t )] 4 



which according to fl38l) is to be evaluated at r such that (x — x(tq)) 2 = 0. In other 
words, in this case signals propagate purely along null intervals: the hyperboloid we 
had for z m > converges to the lightcone, and the retarded time ( 139]) is now at the 
point of intersection between the quark worldline and the past lightcone of x, i.e., 

*ret = t~ \x - f (tret) I- 



4 Some Applications 
4.1 Harmonic motion 

As a first application of our general result ( 13TT) . we will study the example of a 
quark undergoing one-dimensional harmonic motion. This is the setup where, under 
a linearized approximation to the string dynamics, the authors of [5] obtained a 
propagation pattern with a broad tail, seemingly in conflict with the no-broadening 
result of [151 EH] an d the present paper, but, on the other hand, in consonance with 
general expectations for a strongly-coupled non-Abelian system HHJ EES] . In what 
follows, we will show in particular that, despite appearances, the result reported in 
[5] correctly reproduces the linearized limit of our own result. 

Let x = x l denote the direction of oscillation. The trajectory is then given by 

x At = (t,Asm{ut),0,0) , (44) 

with corresponding four-velocity 

«' t = (7,7wAco8(wt),0,0) , 7 = 7= = = . • (45) 

a/1 — bJ l A l COS 2 {UJt) 
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We will consider first the case (not examined by [5j [151 CSj) where the quark 
has a finite mass. According to (13 ip . determining the gluonic profile in this case 
requires, beyond the position, velocity, acceleration and jerk of the quark, knowledge 
of the total external force applied to it, together with its first and second derivatives, 
sometimes called yank and tug. The generalized Lorentz-Dirac equation ffTTj) specifies 
the four-force = j(F ■ v, F) that corresponds to any given quark trajectory. For 
one- dimensional harmonic motion, this equation demands that F = F 1 satisfy 

dF , — -r- — I T~F uj 2 Asin(ut) fl-z±F 2 \ lA , 

— = -Jl-u 2 A 2 cos 2 (out )Jl - z^F 2 - v ' . ^ — . 46 

dt v V z m 1 -u 2 A 2 cos 2 (tut) V z 2 m ) 

Notice that the nonlinearity of this equation implies that the force is in general not 
harmonic, as a consequence both of the extended nature of the quark and of the 
damping effect due to the emitted radiation. 

In the textbook analysis of the forced and damped harmonic oscillator, the forcing 
is prescribed to be harmonic and the oscillator motion is synchronized with this forcing 
only after the decay of a transient component. Conversely, purely harmonic motion 
would be associated with a force that is initially not harmonic. Numerical exploration 
of (I46p reveals similar behavior: if we prescribe the motion of the quark to be purely 
harmonic, as in (I44p . then for generic initial condition, the force contains a transient 
component that dies down in a time interval of order z m . At late times, the force is 
found to be purely oscillatory, but not quite harmonic. Now, as shown in [191 120] and 
recalled in Section [21 the dependence of the string embedding on the external force 
codifies an infinite number of higher derivatives of the quark's trajectory, so the initial 
condition supplied to the numerical integration of (1461) at any finite time should in 
fact be deduced from the behavior of the quark at previous times. If we truly want 
to explore the situation where the quark has been undergoing harmonic motion at all 
times, then we are automatically forced to work in the late-time force regime where 
there is no transient. Equivalently, we must fine-tune the initial condition so that the 
force is purely oscillatory from the start. 

To use equation (jHTj) . we also need to solve (13S]) or its noncovariant equivalent (I3"9"j) 
to obtain the retarded time r or t ret at which all dynamical data are to be evaluated. 
For harmonic quark motion with arbitrary oscillation amplitude, this equation is 
transcendental, and must also be solved numerically. Using the information from both 
numerical integrations in our general expression ( 13T|) . we have obtained the gluonic 
field profile shown in Fig. 4, which shows that the oscillatory motion generates a 
non-linear wave, with crests splitting off from the gluonic cloud of the quark every 
half-cycle. As stipulated by f l3T|) . the overall pattern decays very fast (oc l/|a;| 4 ), 
showing no sign of the characteristic radiation falloff. The waves seen in Fig. 4 are 
thus fluctuations in the near-field of the quark. 

For infinite mass, the quark becomes pointlike and radiation damping becomes 
negligible. The equation of motion (l4"6"j) then simplifies drastically, and implies that 
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Figure 4: Gluonic profile of a heavy quark undergoing harmonic motion in the x-direction, 
with the horizontal axes in units of the quark Compton wavelength, z m = 1, and the vertical 
axis in units of \/A. We choose sample values A = 1 (of the same order as the quark size) and 
u = 1/2. The plots correspond to four successive snapshots of our simulation, for t = 0, 2,4 
and 6, respectively (all within one oscillation period). The finite but proportionately large 
near-field of the quark is capped off in order for the vertical axis to remain at a scale where 
a non-linear wave pattern is clearly visible, with crests splitting off from the gluonic cloud 
of the quark every half-cycle. 

(up to an additive integration constant) 

Fit) = — arctanhL4u; cos(o;t)] , 

i.e., the required external force is purely oscillatory, but (for arbitrary quark oscilla- 
tion amplitude) still not harmonic. Nonetheless, we know that for a pointlike quark 
our formula for the gluonic profile reduces to (145]) . which depends on the assumed 
harmonic motion of the quark, but is independent of the external force. 

As promised in the introduction, we will now show how to make contact between 
our results an those of [5], which likewise determined the expectation value of TrF 2 
for an infinitely massive quark undergoing harmonic motion. The authors of that 
paper carried out their calculation directly in the static gauge a = t, a 1 = z, but of 
course, the fact that we use a different (more geometric) choice of coordinates here 
cannot affect the final result. A difference that does matter is that [5] considered 
only the case where the quark oscillation amplitude is small, A C 1/w, and thereby 
treated the string dynamics in a linearized approximation around the static vertical 
embedding. By inspection of the solution given in j5] , or equivalently, from the full 
solution f|T41) . we see that the fluctuation on the string (together with its derivatives) 
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remains small only for z <C 1/u;, which thus stipulates the region of validity of the 
linearized approach. Through the UV/IR connection, this implies that the result of 
[5] for the gluonic profile is quantitatively reliable only for observation points that 
are close to the source, in the sense that \x\ < Altogether then, we will want 
to examine our general SYM result in the small-amplitude and small-distance (or, 
equivalently, low-frequency) regime 

ujA < 1 , uj\x\ < 1 . (47) 

When working in the linearized regime, it is convenient to rewrite the quark tra- 
jectory (144 p in terms of an exponential, x M = (t, Ae~ luJt , 0, 0), with the understanding 
that henceforth one should in the end take the real part of any complex expressions. 
This notational change (and shift of phase from a sine to a cosine) will also facilitate 
the comparison with the results reported in [5]. Given our general result ( 143]) for 
a pointlike quark (which, it is worth emphasizing, does not depend on the external 
force), the one element we need to obtain an explicit expression for the gluonic field 
is the retarded time defined by the null interval condition t rct = t — \x — x(t Ict )\. To 
linear order in A, it is easy to see that the solution is given by 



^ret 



At 

t - If I + 2± e -™{*-\*\) + o(A 2 ) . (48) 
\x\ 



Using ( I48p in ( I43p . and taking into account the first of the linearization conditions 
( 14T|) . the profile of the gluonic field is found to be 

{TlF2{x)) = le^W + ^W 1 ^"' 1 " 11 ~ MA) + 0(A2) ' m 

where we have defined A = (A, 0,0). As expected, the zeroth-order term in (T49l) 
correctly reproduces the static contribution [6]. The term linear in A represents the 
dynamical part of the field arising from the oscillatory motion of the quark, and, just 
like the full profile (1431) . manifestly contains a single time delay, corresponding to 
propagation strictly at the speed of light. 

Let us now revisit the results of [5]. From the linearized solution for the string 
embedding, the authors of that paper derived the static field contribution shown in 
the first term of (149!) (known previously from [6])@, and going beyond that, deduced 
the leading dynamical contribution to the gluonic field for small oscillation amplitude, 



(TrF 2 (x))g n = r \z'z'\\ - iuz')fWz* + |f|2) e -Mt-v^W-') 



(50) 



with 



r , , 105 57iu 12u 2 iu 3 

IW = — g 7~ + ~ 

u a u a u' u° 



6 It should be noted that Eq. (23) of [5] erroneously reports a denominator that is twice as large 
as that in (l49l). 
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The combination in the integrand of the growing factor z' 2 (l — iuz') with the rapidly 
decaying f(\/z' 2 + \x\ 2 ) ensures that the main contribution to the integral is sharply 
localized around z' ~ \x\, in accord with the UV/IR connection. In (1501) we also see 
explicitly the advertised superposition of contributions with all possible time delays, 
seemingly in conflict with f|49|) . 

To make an honest comparison between expressions fH9|) and (150]) . we would need 
to carry out the integral in the latter. It is important to recall at this point that, even 
though this integral runs over all values of the bulk radial coordinate, the integrand is 
only accurate for z' <C 1/cj, because, as we explained before, the linearized approxi- 
mation breaks down at large values of z' . The result of the integral can therefore only 
be trusted quantitatively at \x\ <C l/oo, so in what follows we will confine ourselves to 
this observation region. Along the lines of [5] , it is convenient to change the variable 
of integration in (1511 to £ = a/1 + z' 2 /\x\ 2 + z' /\x\, to arrive at 



with 



^ F W=i^M6 Z / dtme-**™'" (51) 



MO = jf^& - Mm 2 ~ 1)) 



105(2Q 3 57icu\x\(2Q 2 12a; 2 |x| 2 (2Q . 3 3 " 

L(£ 2 + i) 3 ~ (e + i) 2 (e+T) + lUJ N . 



This form emphasizes the presence of components with all possible propagation ve- 
locities l/£. The function h is a rational function of £ and a polynomial in the small 
quantity u\x\, whose leading term is of order (cu |af| )°. If we integrate by parts, differ- 
entiating the exponential and integrating h(£), we are left with a surface term of this 
same order, which is evaluated at the lower endpoint £ = 1 (with the upper endpoint 
giving a vanishing contribution), and a remaining integral of order (u^o?!) 1 : 



(TrF 2 (x)) dyn 



V\(x ■ A) 

8tt2| 



S7THX 



— {1280 + 15z(19tt - 6A)u\x\ - 1016(a;|f|) 2 

+2^(128 - 57r)(w|x|) 3 + 40(w|f|) 4 } e'^-^ 



— ioj\x\ / d^H(^)e 



-iui(t— £\x\) 



(52) 



where H of course denotes the primitive of h. 

From (|5"2"1) we see that, to order (a;|x|) , the time-dependent gluonic field reduces 
to ^/\{x ■ A) / Att 2 \ x| 6 , in precise agreement with ( 149|) . We can continue integrating 
by parts iteratively to completely reexpress the result of [S] as a series in powers 
of the small parameter uj\x\. Already at order (wlx]) 1 one finds disagreement with 
the exact small-amplitude result ( I49p . consistent with the fact that the integrand in 
( 150]) is not quantitatively reliable away from the small- distance regime (T4T]) . In other 
words, carrying out the worldsheet integral and then linearizing, as we have done in 
this paper, does not yield the exact same answer as linearizing and then integrating, 
as was done in [5]. Nonetheless, what is clear is that, at any order, upon carrying out 
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the integral in (150]) one inevitably evaluates the result at the endpoint £ = 1, and is 
therefore left with a single phase factor exp[— iu(t — \x\)]. 

The final conclusion is thus that the integrated expression (j50j) . which manifestly 
displays the gluonic field as a sum of contributions with all possible time delays, 
actually gives the same result (within its range of validity) as the linearized and 
infinitely-massive version of our formula fl3T|) for the gluonic profile, which incorpo- 
rates a single time delay. This proves that, despite appearances, there is in fact no 
contradiction between [5] and [T^JUS]: it is just that the choice of worldsheet coordi- 
nates and the lack of an exact solution for the string embedding prevented the authors 
of [5] from being able to carry out the integral explicitly to find the net retardation 
pattern. 

4.2 Uniform circular motion 

Let us now apply our general result (|3T]) to a quark undergoing uniform circular 
motion. This is another natural benchmark example, whose associated energy density 
distribution (T 00 (x)) has been studied in detail for the case of an infinitely- massive 
quark [14] (see also [25]). Surprisingly, the radiation field was found to coincide 
with the synchrotron pattern familiar from classical electromagnetism, and therefore 
gave the first indication of the absence of temporal/radial broadening emphasized in 
[15] ITS] . In this section we will consider the more general case of a quark with finite 
mass, and examine the gluonic field via the observable (TrF 2 (x)), which we already 
know to be sensitive only to the near-field of the quark. 

The trajectory of a quark moving with constant angular velocity u along a circle 
of radius p m can be described by 

X 11 = (t, p m COS U)t, p m SU1 Ojt , 0) , (53) 

with corresponding four-velocity 

V ^ L = (7, -l^Pm sin ut, jup m cos cut, 0) , with 7 = = . (54) 

V 1 ~ U Pm 

As in the previous section, to use our finite-mass formula ( 131]) expressed in terms of 
physical quark data via ( 132]) -( 13~T]) . we would need first to solve the Lorentz-Dirac-like 
equation ( |17p to determine the external four-force acting on the quark. 

Fortunately, for uniform circular motion there is a shortcut: it is easy to directly 
deduce the value of the auxiliary tilde variables (describing the motion of a fictitious 
string endpoint at the AdS boundary) that appear in fl3TT) . The reason is that in 
this case the string embedding (115)) or ffT6]) is known to take the form of a uniformly 
rotating spiral [H], and so, to satisfy the requirement f lTT]) that its physical endpoint 
at z m undergo uniform circular motion, we merely need to require that the fictitious 
endpoint at z = also rotate uniformly (with an appropriately shifted phase). Since 
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the radius of the spiral grows with the radial AdS coordinate according to [T] 



P( z ) =Po\h- + 



U 2 Z 2 



1 - u 2 p 2 ' 

we must choose the radius at the boundary, p Q , such that 



POm) = p \ 1 + 



u 2 z 2 



u) 2 pl 



Pr, 



(55) 
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Figure 5: Gluonic profile of a heavy quark undergoing uniform circular motion in the xy- 
plane, with parameters z m = 1, p m = 1 and oj = 1/2 (for which (i55j) yields /9o=0.874). 
The plots correspond to four successive snapshots of our simulation for t = 0, 2, 4 and 6 
respectively (all within one period). The near-field waves display no broadening, and are 
clearly consistent with the radiation pattern observed in \14:\ [26] [To] . 

To derive explicit results we additionally need to determine the retarded time 
tret by solving ( )39l) . Just like in the previous section, the equation turns out to be 
transcendental, so we process it numerically. The final result for (Tr F 2 (x)) is shown in 
Fig. 5, for a sample choice of parameter values. The spiral profile of the gluonic fields 
is clearly in close resemblance with the pattern detected in [T4"| 126] [15] through the 
observable (T 00 (x)). Indeed, one can check that as the quark of the mass approaches 
infinity, the pattern precisely reduces to the expected synchrotron form. At finite 
mass, we have a generalized synchrotron pattern, with somewhat different width of 
the spiral arm and a net subluminal propagation speed, as dictated by (|39l) . Still, no 
temporal broadening is observed and the disturbance remains sharply localizedjj 

7 It would be interesting to follow the ideas presented for the energy density distribution in [26] to 
try to recover these same results in terms of a superposition of dilatonic shock waves, but we leave 
this for future work [34] . 
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4.3 Constant acceleration 



As a final example, let us consider a quark moving with constant proper acceleration. 
The nature of the radiation emitted by a uniformly accelerated charge in the frame- 
work of classical electrodynamics is a topic that has received considerable attention 
over the years (see for instance j35j and the references therein). Various authors have 
already used holographic techniques to study this configuration in the previously un- 
charted, strongly-coupled regime [221 ESI [21] • In particular, it was shown that the 
physics behind the Unruh effect (on the field theory side) is nicely captured by the 
structure of the string embedding from the gravitational viewpoint^] The purpose 
of this subsection is to complement the previous analyses by showing explicitly the 
gluonic pattern sourced by such a quark. 

Without loss of generality, we will assume that the motion is along the x-direction 
so, in covariant notation, the trajectory we are interested in is given by the usual 
hyperbola 

x" = (t, A~Wl + AH 2 , 0, 0) , (56) 
where A is the magnitude of the four acceleration, 

As \l¥¥=£™- (57) 

and where for convenience we have made a particular choice of the spacetime origin. 
From f )56|) it follows that 

7 (t) = Vl + AH 2 , (58) 

and 

r = A _1 arcsinh(At) , (59) 

where r denotes, as usual, the time measured by the internal clock of the accelerated 
observer. The worldline described by such a hyperbola is displayed in Fig. |6j 

Again, determining the gluonic profile in the general case (1311) requires knowledge 
of not only the quark trajectory and its derivatives, but also the total force applied 
to it. In this case, the equation of motion (fTTI) reduces simply to 

A(l - z^F 2 ) = ZmF^l-z^F 2 + z 2 m FVl + AH 2 . (60) 

Notice that, just as in the examples that we have studied before, the equation obtained 
is a non-linear first order differential equation, and this automatically implies that 
there are an infinite number of solutions depending on the initial condition. Numerical 
integration of (160)) reveals a structure similar the one found in Section 14.11 for generic 
initial conditions, the solution is a non-linear superposition of a transient response 
that decays rapidly (within a time scale of order z m ) and a forced (steady-state) 
solution that dominates at late times and is independent of the initial condition. 

8 The same mechanism was generalized recently to the case of an arbitrary accelerated trajectory 
on a bounded space [36] . 
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Figure 6: The hyperbolic trajectory x^(t) of a quark under uniformly accelerated motion. 
The retarded time t re t associated with any point is determined by the intersection between 
the quark worldline and the hyperbola (138j) contained inside the past light-cone of the given 
point. Only regions I and II are affected by the accelerated quark. 



This is presumably an inherent property of the structure of the generalized Lorentz- 
Dirac equation (TTTj) . associated with the presence of radiation damping, and, given 
the acquired numerical evidence, we expect it to hold also for more general/arbitrary 
trajectories. 

As noted in Section I4.1[ the choice of initial condition for the force encodes the 
past behavior of the quark, so if we wish to describe the situation where the quark 
has been undergoing uniform acceleration for all times, we must effectively ignore the 
transient tail. Equivalently, we can give the initial condition for the force at a finite 
time, but tune it to coincide with the asymptotic, steady-state, value. For arbitrary 
A, this value is seen from ([60]) to be 



If our interest lies solely on the perpetually accelerating quark, we could have 
in fact avoided any mention of the Lorentz-Dirac-like equation (160]) . because in this 
case the string embedding (|T5| or (|T6|) is known to take the simple form X(t,z) = 
■s/A~ 2 + t 2 -z 2 + zl [22] [23], which makes it easy to directly deduce the behavior 
of the auxiliary tilde variables required in (13T~|) : uniform acceleration at the physical 
string endpoint z = z m corresponds to uniform acceleration (of different magnitude) at 
the fictitious endpoint at z = 0. We should stress that this simple relation between the 
behavior of the string embedding at different radial depths, which we also encountered 
for the case of uniform circular motion in Section 14.2] is not by any means generic, 
and in particular does not hold for the oscillating quark studied in Section 14.11 (for 
which harmonic oscillation of the physical endpoint is easily seen to translate into 
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nonharmonic oscillation of the fictitious boundary endpoint). 

To determine the gluonic profile through (1311) . we additionally need to solve (I3"8"j) 
to obtain the retarded time r , and then evaluate all the data at this time. For con- 
venience we will use its noncovariant equivalent fl39|) . which can be inverted explicitly 
to deduce that 



This can be interpreted geometrically as follows. The retarded time t ret associated 
with a given point (t, x, y, z) corresponds to the (unique) intersection of the hyperbola 
fl3"9j) . contained inside the past light-cone of (t, x, y, z), with the quark worldline x^{t). 
For instance, in Fig. [SJ the points p' and q' define the retarded times t^ et and t^et 
associated with the points p and q, respectively. The fact that the signal travels along 
the hyperbola ( |39l) instead of along a null line means that the velocity of propagation 
is subluminal; we will come back to this point in the next section. For z m — > the 
hyperbola converges to the lightcone, and the retarded time is determined by the 
point of intersection between the quark worldline and the past lightcone. In this case 
the signal travels strictly at the speed of light. Note also that only regions I and II 
are affected by the fields due to a charged particle with a worldline given by ([56]) . On 
the contrary, points in regions III and IV, are causally disconnected from the quark 
worldline (see for example the point s in Fig. E]) and this automatically implies that 
the gluonic field vanishes there. 

Finally, using (|56|) and (!6T|) and evaluating (I3~TI) in (1621) we obtain the gluonic field 
profile shown in Fig. [7J We choose sample values z m = 1, \f\ = 1, A = 1/2 and 
we plot points in (part of) regions I and II. The figures correspond to four successive 
snapshots of the profile, for t — 2, 4, 6 and 8 respectively. The most salient difference 
of the finite mass case with respect to the z m — > case is that the extended nature 
of the quark makes the gluonic profile finite even at the location of the quark. The 
overall pattern decays very fast however, and for points very far from the quark there 
is no significant difference between these two cases. The behavior of the gluonic field 
displays the expected Lorentz contraction in the longitudinal direction as the quark 
increases its velocity. 

It would be interesting to investigate this setup further through the calculation 
of (Tjj^x)) for both the inertial observer and the comoving observer, but this com- 
putation lies beyond the scope of this paper. The radiation field (i.e., the part that 
decays as l/|x| 2 ) could be extracted from that observable, and contrasted with the 
expected result from classical electrodynamics [53] (with which contact has already 
been made for the energy density in the case of an infinitely massive quark [IS]). It 
would also be interesting to exploit the tensorial properties of our result as in J37J [23] 
to obtain the gluonic profile as seen by the comoving observer. 




(62) 
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Figure 7: Gluonic profile of a heavy quark with uniform proper acceleration in the re- 
direction, with the horizontal axes in units of the cjiicirk Compton wavelength, z m = 1, and 
the vertical axis in units of vA- The plots correspond to four successive snapshots of our 
simulation, for t = 2,4,6 and 8, respectively, and acceleration A = 1/2. The value of the 
gluonic profile at the location of the quark is finite but this section is capped off to show 
the near field details. The pattern obtained displays the expected Lorentz contraction in 
the transversal direction as the quark increases its velocity. 

5 Discussion 

Let us now discuss the physical implications of our results. According to ( I3TT) . data of 
the quark (position, velocity, acceleration and jerk, together with external force, yank 
and tug) at the sole instant fl39l) on its trajectory determine the expectation value of 
Tt F 2 (x). It is important to remember that, just as in the equation of motion (fTT|) . 
the presence of the external force in our final expression for the one-point function 
conveniently replaces a dependence on an infinite number of higher-derivative terms 
[20] . and is therefore indicative of nonlocality. In our result, these higher derivatives 
appear modulated by powers of the quark Compton wavelength z m , so, in essence, 
it is data within z m of t rot that have an impact on the gluonic field at the chosen 
observation point. The nonlocality is thus naturally associated with the extended 
nature of the quark. 

This nonlocality disappears in the pointlike/infinitely-massive limit z m — > 0, 
where, as noted at the end of Section [3], our general result reduces to fH3|) . which 
depends only on the position and velocity of the quark at a retarded time evaluated 
on the past lightcone of the observation point, just as in classical electrodynamics. 
As a matter of fact, our result ( ]43i) for an infinitely massive quark can easily be 
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seen to agree with the square of the Lienard-Wiechert electromagnetic field strength 
(more precisely, (1/4)F 2 ) set up classically by a pointlike electric charge e [TJ, under 
the replacement^ e 2 — > \/~\/8ir 2 . Equivalently, it agrees with the leading order result 
within M = 4 SYM at weak coupling, with the rescaling A — > 4a/A, which is exactly 
the same connection found in [16] for the expectation value of the energy density. 
As in electromagnetism, the gluonic field ( l4"3]) is simply the boosted version of the 
Coulombic profile set up by a static pointlike quark [5], at the appropriate retarded 
time. 

In the pointlike (z m = 0) case, then, the retardation pattern found for the gluonic 
field observable under present examination, (Tr F 2 (x)), is seen to be exactly the same 
as that obtained for the energy density, (T 00 (x)), in [H] and, especially, [151 EE] ( as 
well as for other observables in pj)]): the net profile propagates strictly at the speed 
of light, with no temporal (or, equivalently, radial) broadening. 

In the case of a non-pointlike/finitely-massive (z m > 0) quark (which was not 
examined in [HI [151 EB]), we have found that two features of the gluonic profile are 
modified. First, signals no longer travel at the speed of light, but rather propagate 
along the timelike interval ( 13"H|) . The average velocity of propagation of the signal, 

^average = (x - X (tret))/ ~ *ret) follows from §M) as 



z 2 



^ average — -\ / 1 [t — t ) 2 ' 

Since z m < t — t TCt < oo, (1631) implies that < leverage < 1- From the gravity perspec- 
tive, (1391) clearly just means that the net signal from the string to the observation 
point on the AdS boundary has the smallest possible delay: the time for null propaga- 
tion starting at the string endpoint (as opposed to any other location on the string). 
The interpretation of this delay from the gauge theory perspective is more interesting: 
it is due to the finite size of the dressed quark. Consider the seemingly most peculiar 
case where ^average = 0, which happens when t — t Tet = z m . In this case x — x{t rct ) = 0, 
so the spatial location of the quark at the relevant instant coincides with that of the 
observation event, and it might seem surprising that there is still a delay before the 
gluonic field 'takes notice'. Recalling, however, that z m is nothing but the size of the 
quark, we see that this result makes perfect physical sense: the delay just indicates 
the time for information to propagate from the 'edge' of the quark to the observation 
point located at its (instantaneous) center. We can see in ( 13"§j) that this same delay 
is present in the more general case where leverage > 0, but becomes progressively less 
noticeable as the distance between the quark and the observation point increases. 

We emphasize that the full pattern of propagation is very complicated, with each 
point on the non-Abelian medium potentially yielding a contribution to the gluonic 
field at the observation point, and t> av erage is associated solely with the delay in the net 
result. At first sight, it might seem odd that this velocity is not fixed, but this is in 



9 Notice that this is somewhat different from the replacement e 2 — > 3V~\/4:Tr that translates the 
radiated energy in classical electromagnetism to the one in strongly-coupled SYM [17] . 
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fact unavoidable: for subluminal propagation, the only Lorentz-covariant possibility 
is for signals to be emitted along hyperboloids at constant timelike interval, precisely 
as prescribed by (|3"H]1 . 

A second modification in the non-pointlike profile ( 131]) with respect to (143)) is that, 
as we mentioned at the beginning of this section, it depends on more data than just 
the position and velocity of the quark. In particular, the acceleration, jerk and snap 
of the quark have an impact on the gluonic field. Given the spatial falloff (oc l/|x| 4 ) of 
(TrF 2 (x)), these higher-derivative contributions cannot be interpreted as indicative 
of the presence of radiation. Of course, radiation is expected to be generally present 
when the quark accelerates, but its unambiguous detection requires a generalization 
of the calculation in (Too(x)) to the case of finite mass. The higher-derivative terms 
in our result for (TrF 2 (x)) are due to (or rather, evidence) the fact that the gluonic 
cloud of the quark is deformed when the latter accelerates (or jerks or snaps) |2Uj . 

In spite of these differences, one important feature shared by the pointlike and non- 
pointlike results is that they are determined by the behavior of the quark/endpoint 
at a single retarded time tq or t ret . It is worth recapping here how this result came 
about. At the start of the computation, we had the expression (fT2"j) for the dilaton 
field in the bulk, which is a priori sourced by all points on the string worldsheet. As 
we emphasized below ( 124)) and illustrated in Figs. 2 and 3, the causal structure of 
the propagator ([8]) appropriately restricts the source points to lie within a specific 
swath on the worldsheet, which covers a finite temporal extent (parametrized by t', 
t' r or r') at any given radial depth in the bulk (parametrized by z' or U). As we 
also noted, this swath shrinks to a line when, in accordance with the GKPW recipe 
[3], the observation point is taken towards the AdS boundary. The requirement of 
propagation all the way to z — > thus picks out only one instant on the worldsheet 
at each radial depth, and the dilaton/gluonic field profile at a given (Minkowski) 
spacetime location is obtained as the integral over all such contributions. The final, 
crucial ingredient in our calculation was the fact that, as seen in (l2"3"j) . the integrand 
in this radial integral, when expressed in terms of the (r, U) parametrization, turned 
out to be a total derivative, allowing the result to be written purely as a surface 
term. This is why, when all the dust settles, the observed (TrF 2 (x)) depends only 
on the behavior of the (lower) string endpoint (with the upper endpoint excluded by 
causality), and thus involves a single retarded time. 

In the calculation, our being left with a total derivative in (1231) can be traced back 
to the cancelation of the z' 2 terms in (|2~T1) . An analogous cancelation was seen to take 
place in the computation of the energy density [16] . Its appearance here in the midst 
of a scalar calculation shows that the effect is not specific to the interplay between 
the various tensor components relevant to [TB]. In our setting, the cancelation clearly 
stems from the specific structure of the string embedding (TT41) or (TTBT) . Physically, 
then, the absence of temporal broadening in our result is closely associated with 
the use of a purely retarded solution. The same feature is seen when using a purely 
advanced configuration, but would not be expected in the more general case involving 
nonlinear superpositions of retarded and advanced solutions (or in the presence of 
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a thermal medium [12], which, after an initial period, naturally brings into play 
analogous superpositions [TBI I3"3"] ). 

Given that we are dealing with a strongly-coupled non-Abelian gauge theory, it 
is remarkable that we have obtained a gluonic field profile that can be expressed in 
the relatively simple form (|3T|) (or, in the pointlike quark limit, (I43p ). which depends 
on a single retarded time ( 138]) . Since the Yang-Mills field is a nonlinear medium, 
we would naturally expect it to generate infinitely broadened disturbances, arising 
from rescattering at all possible length/energy scales of the original signal emitted 
by the quark, leading to a pattern such as the one reported in [5], with components 
propagating at speeds arbitrarily slower than light. Indeed, the authors of [151 US] 
have argued that the absence of temporal broadening in results such as (H3|) is a 
physically unwanted feature, and indicates a failure of the supergravity approximation 
of AdS/CFT to capture the full quantum dynamics. 

In our computation, however, we have seen how the final result ( }3~T|) or fj43l is 
in fact assembled from contributions of all points along the string, which, via the 
UV/IR connection [7], translate precisely into the expected contributions from the 
SYM fields at all possible length scales. The appearance of a total derivative in (|2"3"|) 
enabled us to explicitly express the result of the radial integral as a boundary term, 
but this does not imply that all the contributions from interior points of the string are 
canceling out, only that their sum can be rewritten purely as a (different) function of 
data at the string endpoint. This property would have been much more difficult to 
notice if we had chosen to carry out the analysis in the original static gauge (t',z'), 
or, more generally, in a parametrization not closely related our geometric choice of 
the retarded proper time r' and the invariant AdS distance U . In such cases, we 
would be left contemplating a sum over contributions that manifestly involve time 
delays corresponding to all possible subluminal speeds, even though, of course, the 
final result of the sum is reparametrization invariant and would therefore necessarily 
coincide with what we have found here. As we have seen very explicitly in Section 
14. 1} the study of an oscillating quark in [5] provides us with a pre-existing example 
of precisely this sort of situation. 

Our physical interpretation is thus different from that of [151 ES] > even though the 
properties of our results are in full agreement with theirs. We see the AdS calculation 
yielding a one-point function that is manifestly obtained as a sum over contributions of 
the gluonic field at all possible length scales, just as one would expect, and therefore 
take f l3T|) and f l4^|) at face value, as a prediction of the AdS/CFT correspondence 
for the spacetime pattern emerging in the A, N c — > oo limit where the supergravity 
calculation is performed, for the specific case where the field configuration is assumed 
to be purely outgoing. As we mentioned in the Introduction, at least in retrospect it 
feels to some extent natural for this assumption to imply a restriction on the overall 
retardation pattern. What is certainly surprising is the actual form of the predicted 
restriction: even though we would a priori expect the net result to depend on data 
characterizing the entire quark trajectory, our final pointlike profile ( 14U]) is found to be 
completely local, and its non-pointlike extension (|3T|) is seen to have a nonlocality that 
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(being bounded by the quark Compton wavelength (JSJ) is directly associated with the 
extended nature of the finitely-massive quark, rather than with the nonlinear nature 
of the SYM fields. 

There are several directions along which this circle of ideas might be profitably 
extended. It would certainly be interesting to inquire to what extent the surprising 
pattern of unbroadened propagation obtained in [15], [16] and the present paper is 
present in other setups made available to us by gauge/gravity duality (e.g., in non- 
conformal field theories). Another such direction, which we hope to report on in the 
near future [38], is to revisit the calculation in [16] of the energy density in the gluonic 
field generated by a quark undergoing arbitrary motion, and generalize it to the case 
of finite quark mass. Taken together with the results of the present paper, this could 
shed light on the appropriate spacetime separation of the near- and radiation- field 
contributions, examined previously from a worldsheet perspective in [T8"| Q3JJ [2D] , and 
also touched upon in [16]. Through these and a myriad other ongoing explorations, 
the correspondence still seems poised to teach us many new and important lessons 
about the dynamics of strongly-coupled gauge theories. 
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A Appendix 

Here we give some more detail on the calculation of the one-point function of Tr F 2 
for an arbitrary quark trajectory presented in Section [3j It helps to summarize the 
results (TJZD, (J29J, §Q by writing 



where the coefficients Aj are independent of U m \ D and z, and are given explicitly by 



f = f + A 1 U min z + (A 2 + A 3 UlJz 2 + 0(z 3 ) , 



(64) 
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z m + (x- x') ■ (v' + a!z m ) 
1 



A 2 



2 [z m + (x- x') ■ (v 1 + a'z m )\ ' 
z 2 m (l + (x-x')-(a' + ~fz m )) 
2 [z m + (x - x') ■ (v' + a'z m )f 



(65) 
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Straightforward algebra then shows that 

(x - x{t')) 2 + z 2 = B 1 + B 2 U min z + (B 3 + B 4 U 2 lin )z 2 + 0(z 3 ) , (66) 

with 

B 1 = (x- x') 2 , 

B 2 = -2Ax(x - x') ■ v' , (67) 

B 3 = 1 - 2A 2 {x - x') ■ v' , 

Bi = -2A 3 {x-x')-v'-Al{l + (x-x')-a') ; 

whereas 

^ m + (x - x{f')) ■ (v(f) + z m a(f')) = d + C 2 C/ min z + (C 3 + C 4 C/ min )^ + 0(2 3 ) , (68) 
with 

C\ = z m + (x - x) ■ (v 1 + az m ) , 

C 2 = A 1 (l + (x-x')-(a , + ~j'z m )) , (69) 
C 3 = A 2 (l + (x-x')-(a'+~j'z m )) , 

C A = A 3 (1 + (x - x') ■ (a' + ~fz m )) + ~A\ ((x - x') ■ Q' + s'z m ) + z m a' 2 ) . 

Combining ( 12"B"|) . ( 16"6"j) and (I6"8"j) and expanding df' z/((x — x{f')) 2 + z 2 ), we find 
terms of order z 2 and z 3 that have up to a linear dependence on ?7 mm , and therefore 
vanish upon carrying out the integral fl25l) . The leading contribution to the dilaton 
field in the z — > limit arises from the terms of order z 4 ?7 mm , which yield 

_ V\z m z 4 ( B\{C 2 - CiC 4 ) + B X B 2 C X C 2 + (£?j - £?i£? 4 )C 2 ^ 5 , 
32tt 2 V BfCf )+ u V)- 

(70) 

Plugging this into ([7]) and using fl67|) and (169]) . we finally arrive at the result fl3Tj) . 
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